Geometry and elasticity of strips and flowers 



O 
O 
(N 

Oh' 

(D 



o 



o 



> 
in 
m 

On 
O 

O 



O 

o 



M. Marder 

Center for Nonlinear Dynamics and Department of Physics The University of Texas at Austin, Austin TX 78712, f/S/Q 

N. Papanicolaou 

Department of Physics, University of Crete, and Research Center of Crete, Heraklion, Greec^ 

We solve several problems that involve imposing metrics on surfaces. The problem of a strip with a linear 
metric gradient is formulated in terms of a Lagrangean similar to those used for spin systems. We are able to 
show that the low energy state of long strips is a twisted helical state like a telephone cord. We then extend 
the techniques used in this solution to two-dimensional sheets with more general metrics. We find evolution 
equations and show that when they are not singular, a surface is determined by knowledge of its metric, and the 
shape of the surface along one line. Finally, we provide numerical evidence that once these evolution equations 
become singular, either the surface is not differentiable, or else the metric deviates from the target metric as a 
result of minimization of a suitable energy functional . 
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I. INTRODUCTION 

In a series of experiments performed in 2002, Sharon et. 
a/.fl] showed that thin sheets deformed in a smooth fash- 
ion spontaneously fold into convoluted shapes with much less 
symmetry than the original deformation. These experiments 
raised many questions about the relationship between local 
changes in distance within a sheet, and the global shape the 
sheet adopts. The first theoretical papers that set about ex- 
plaining the experimentsJ3 |3l [sl pl established certain 
points, but left many others unsettled. What was most clearly 
established was that one should study the problem by choos- 
ing a metric function suggested by the physical process that 
deformed the sheet, and trying to determine what surfaces 
were compatible with this metric. Much of the effort in these 
papers was devoted to exploring the problem of a long thin 
strip of material with a linear metric gradient in one direction. 
This strip problem was not solved in great generality, and its 
relation to the more general two-dimensional problems was 
not definite. Furthermore, the basic idea of the calculations 
was not clear. On the one hand, the papers spoke of imposing 
a metric on a sheet, while on the other there were many indi- 
cations that the actual metric of the final shape adopted by the 
sheet was something else. 

This paper has three main sections which address the ques- 
tions we have just described. 

First, we return to the problem of the long strip with a met- 
ric gradient and provide a more complete solution than was 
afforded previously. In particular, we find the lowest-energy 
shape of this strip when it becomes very long. In this limit, 
the strip coils into a helix like a telephone cord. The edge of 
a flower or a leaf cannot behave in this way. Therefore, solu- 
tions of the strip problem have only limited relevance to the 
more general two-dimensional questions raised by the experi- 
ments. However, the formal techniques used to solve the strip 



prove to be helpful, and provide the basis for our approach to 
the two-dimensional problem. 

Second, we return to the problem of two-dimensional 
sheets, and attempt to determine the extent to which speci- 
fying a metric determines the surface. We can discard the 
possibility that a metric alone can dictate the shape of a sur- 
face. For example, a piece of paper can be bent into infinitely 
many smooth shapes, all of which share the same Euclidean 
metric. This property is not simply due to the fact that the 
piece of paper is initially flat. Tape the paper into the shape 
of a cylinder, and still it can be deformed in infinitely many 
ways without change of metric. A classic theorem of differ- 
ential geometry states that surfaces are uniquely determined, 
up to rotations and translations, by the metric, and by the sec- 
ond fundamental form, provided that they satisfy the Gauss- 
Codazzi equations! 7]. However, the second fundamental form 
provides more information than is needed. We will show that 
if one knows the metric of a surface, and the shape of the sur- 
face along one line, then the rest of the surface is determined 
by an evolution equation. This statement must be qualified. 
The procedure that determines the surface typically contains 
singularities, and the attempt to find the surface with geometry 
alone comes to a halt. 

Third, we discuss numerical procedures to find surfaces 
corresponding to metrics for which the evolution equations 
of the previous section become singular. We describe in some 
detail an energy functional whose minima should provide sur- 
faces that approach a desired target metric as closely as possi- 
ble. In cases where evolution equations dictate surfaces based 
upon geometry alone, the energy functional recovers them. 
When the evolution equations fail to find a surface, then ei- 
ther the energy functional finds a non-differentiable surface, 
or else it finds a smooth flower-like surface whose metric is 
different from the one that was supposedly imposed. 



n. GROUND STATE FOR TWISTED STRIPS 
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Reference |8] shows a number of thin strips cut from the 
edge of a leaf. Each of them curls up into a circle, with the 
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curvature of each strip depending upon the gradient of the 
metric at that point in the leaf. Such observations suggested 
that the rippHng pattern at the edge of a leaf could be under- 
stood by focusing upon a thin strip of material with a metric 
gradient. Such a strip, freely allowed to seek out its lowest 
energy state simply curls up into a circle. Therefore, the addi- 
tional constraint was added that the two ends of the strip had 
to be some distance A apart|2]. This problem yielded non- 
trivial solutions^ Isl IslQl- Certain special cases could be 
found analytically. Somewhat more general solutions could be 
found numerically. However, the numerical procedures were 
not very stable, and left open a number of questions. For ex- 
ample, there was some speculation that as the length L of the 
curved strip became infinite with L/X kept finite, the lowest- 
energy state might be fractal. This speculation was incorrect, 
as we now show through a correct solution of the problem. 

A. Problem setting 

We briefly review the specification of the strip problem. 
Consider a circular strip of paper, as shown in Figure^ of ra- 
dius R, length L = 2t:R, width w, and thickness t. The length 
L can be arbitrary, and may be greater or less than 2ttR. We 
will be interested in strips where the following limits apply: 

/u;\2 t w L 

[r) 

When these conditions hold, the energy of the strip takes a 
very simple form as a functional of a line passing through the 
center of the strip|2]. 

In particular, let s be the arc length along the center of the 
strip. The low-energy conformations of the strip are captured 
by two orthonormal vectors fi(s) and ^2(5), where fi points 
along the arc s, and f2 lies in the plane of the strip and is 
orthogonal to fi . Define also = fi x 7^2 , so as to obtain 
a set of unit vectors that describes in a natural way the local 
orientation of the strip. The location r(s) of each point along 
the center-line of the strip can be obtained from 

f(s)= / ds'his'). (2) 
Jo 

At each point s' along the way, the precise orientation of the 
strip is specified by the unit vector f2, or equivalently by the 
unit vector fa, which is normal to the strip. Therefore, the 
vectors fi(s) and f2(s) contain all the information necessary 
to deduce the shape of the strip when it is embedded in three- 
dimensional space. The inequalities Q ensure that the strip 
remains essentially flat and featureless in the direction point- 
ing along f2. 

In terms of the vectors fi{s) and f2(s) the energy of the 
strip is given byt^i3|] 

U = [ ds£{s), (3) 
Jo 

where 

S = ^{\r[\'-l/R') + ^m'-l/R'), (4) 




FIG. 1: The buckled strips studied in this section can be examined 
by cutting out this circular strip, of radius R and width w. Roughly 
speaking, the main question addressed here is to predict the shape 
of the strip when the two ends are pulled apart to some specified 
distance and relative orientation. For the general solutions studied in 
this section, the total length L of the strip need not equal 2nR. 



and primes indicate derivatives with respect to arc length s. 

This energy is subject to two constraints. The first is a local 
constraint that results from the fact that the strip has a metric 
gradient in the direction of 7-2 that causes it to want to curl up 
with radius of curvature R. This constraint is captured by the 
condition 

f[ ■ 7-2 = -f'2 ■ fi = l/R. (5) 

The second constraint is that the ends of the strip be pulled 
apart by some distance; this constraint prevents the strip from 
simply curling up into a circle. We impose this constraint 
through 

(f(L) - r(0)) • z = A. (6) 

That is, the z component of the difference between starting 
and ending points of the strip is constrained to have value 
A. By scanning through values of A, one equivalently scans 
through all allowed distances between the starting and ending 
point of the strip. The constraint in the form of Eq. (|6|l is much 
easier to work with formally than if it were literally expressed 
in terms of end-to-end distance. For definiteness, we choose 
a fixed (laboratory) frame as shown in Figure ^ The z axis 
is taken to be perpendicular to the plane defined by the ini- 
tial (undeformed) circular strip, while the x and y axes point 
along the initial directions of fi and f2 at one of the endpoints 
of the strip. 
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The trihedral ri , f2 , fa uses nine coordinates to represent the orientation of the strip, when in fact only three are needed. 
I 



Further analysis is greatly simplified by rewriting the unit vectors in terms of Euler angles (Ref.j^l Eq.4.46), through 



f 1 = [ cos'(/'COS(/) — cos0sin(/)sin^/', cos'(/'sin</) + cos^cosc^sin-i/), sin sin t/']; 
f2 = [— siuT/i COS0 — cos6'sin^cos^, — sin^/isin^ + cos^cos^cosi/', sin 6' cos ■(/'] ; 
fa = [ sin 9 sin — sin 6 cos 0, cos 0] . 



(7) 



r 



With this representation, we can rewrite the local constraint 
Eq. Q as 

cosO <i)' + = l/R- (8) 

again, primes refer to derivatives with respect to the arc length 
s. 

Writing out Eq. in terms of the Euler angles one has 

£ = i [(C2-Ci)cOs2l/. + Ci] 

+ (C2-Ci)cosV'sinV'sin6'0'6'' (9) 

+ i [{Ci - C2) cos2 V + C2] sin^ 6 

In order to impose the global constraint Eq. (|6}, we will em- 
ploy a Lagrange multiplier, and write down the Lagrangean 



[(C2-Ci)cos2v^ + Ci]0'2 
+ (C2 - Ci ) cos ?/' sin iisinO <p'9' 
+ i [{Ci~ 02)005^^ + 02] sin^i 
—hsuiij^smO. 



(10) 



Note that 0' appears in Eq. (I1Q> . but not itself. This fact will 
allow for further simplification. 

The constants Ci and O2 depend upon elastic properties 
of the strip. For a particular case studied in Ref.|5|, they are 
related through 2C2 = 3Ci. Materials with different elastic 
properties would give different constants, so we take Ci and 
O2 as free variables. The structure of the equations is particu- 
larly simple when Ci = C2 , so we will begin with that case, 
and return later to the more general situation. 

Adopting Oi = O2 ~ 1, the Lagrangean of Eq. sim- 
plifies to 



C 



(11) 



2 2 

All appearances of can now be eliminated by employing the 
local constraint Eq. (|8} to give 
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C ^ -9'^ + 2*^^' ^ 1/ RY i&n^ 9 -hsm^ sine. (12) 

Similar Lagrangeans appear in the study of spin svstems full . 
The equations of motion following from Eq. (I12t are 



{ip' - l/R) i&T?9 



9" = (V^' - l/R)' 



sm 



cos^ 9 



= —hcosilJsm9 (13a) 
hsmipcos9 (13b) 



B. Particular solution 

Ref. f?) provided a family of exact solutions, and we search 
for these again with this new formalism. We find that they 
emerge if one sets /i = in Eqs. ( I13> . In this case, one imme- 
diately integrates Eq. ( I13at to obtain 

{ij/ -1/R)tan^9 = P, (14) 

where /3 is an integration constant. Inserting this relation into 
Eq. ( I13bt and integrating gives 



'■~l3^/sin^9 
cos6' = -^1 - (/3/a)2 sin(as), (15) 

where a is an additional integration constant. Without loss of 
generality one can choose a > 0. Note from Eq. (I15> that 
a > (3 if the solutions are to remain real. 
For ip one obtains 

ijj = tan-i ((/?/«) tan(as)) - (/3 - l/R)s + ^po- (16) 
Similarly, for 4> one obtains 



4> ~ ~ tan ^ ^■\/ {a/f3)^ — 1 cos(q:s)^ + 0o; 



(17) 



the integration constants -00 and 0o give the value of i/j and 
when s = 0. 

The global constraint Eq. (|6j requires that 

f-L 

ds sin6'sin'i/' = A. (18) 



As A becomes large, it is only possible to satisfy Eq. (I18> if 
sin 6* and sin?/' have the same period. Comparing Eqs. ( I16> 
and il5i . one sees that this condition can be satisfied by taking 



/3 = a- 
With this choice, one has 



l/R. 



(19) 



sin 9 sin ip — fi 



{aR + sin^ (as)) /{aR), (20) 



which with a bit of manipulation upon setting R = 1 repro- 
duces Eq. (41a) in Ref. |2]; note that the x axis in that ref- 
erence corresponds to the z axis here. The remainder of the 
special solution obtained previously can be recovered as well. 
Since this special solution depends upon setting h — 0, ot 
equivalently to fixing a relationship between A and L not re- 
quired by the original problem, we move to a numerical ap- 
proach capable of solving the problem in greater generality. 



C. Numerical solutions 

We can obtain reliable numerical ground states of Eq. (Ill> 
by changing the form of the Lagrange multiplier to enforce the 
global constraint in a fashion that involves a positive definite 
functional: 



F- 



ds 



H I \ - 



ds sin -0 sin ( 



(21) 



where H is chosen on the order of 100. To find solutions, we 
simply minimize F. 

In the minimization procedure we fix boundary conditions 
on 6 and ip at s = and s = L. We may also choose the 
initial (s = 0) value of but its final (s = L) value is left 
free to adjust according to Eq. (|8}. Now we keep one of the 
endpoints of the strip fixed and orient corresponding initial 
directions of fi and f2 along the x and y axes of the laboratory 
frame. Noting that when 9 = 0, the orientation of the trihedral 
depends only upon + the initial values of the Euler angles 
are set at 



9{0) = 0, 1/^(0) + 0(0) = 0. 



(22) 



At the other endpoint, the values of 9 and ip are chosen arbi- 
trarily; e.g., 



9{L) = 7T/3, ij{L) = 7r/4, 



(23) 



while the end-to-end distance is controlled by the parameter 
A. A numerical solution for L = 13 and A = 0.75i yields 
values for 9 and -ip shown in Figure 121 The angle (p = (j){s) is 




Arclength s 



FIG. 2: Plot of Euler angles 6 and for system with total length 
L = 13, constrained to have height in z direction A = .751/. Note 
that apart from some variation near the endpoints needed to obey the 
boundary conditions, both angles go to constant values throughout 
most of the length of the sample. 



FIG. 3: (a) Three-dimensional visualization of solution shown in 
Figure |2lin color), (b) Visualization of low-energy solution for a 
long strip with \/L = 0.625. The strip winds in helical fashion 
around the z axis (in color). 



obtained by an elementary integration of Eq. (jS), and the com- 
plete orthonormal trihedral from Eq. 0. It is then straightfor- 
ward to calculate r = r(s) from Eq. (|2} which determines 
the relative position of the strip in three-dimensional space, as 
illustrated in Figure|3(a). 

One property of the solution that might seem unsatisfactory 
is that the strip wraps a number of times around the z axis. If 
one were to grab two ends of a strip and pull them apart, this 
would not be allowed. If one wants to find energy-minimizing 
solutions without any net twist around the vertical axis, the 
following strategy is effective: Take a solution obeying the 
boundary conditions 



9{L) =Tr/2, ^(i)=7r/2, 



(24) 



produce a mirror image with z — z, x x, y ^ y. The 
resulting function also minimizes the functional F. It can be 
glued on to the solution found so far, joining smoothly to it at 
s = L because of Eq. MAI , and will continue on to terminate 
at (0, 0, 2A) when s = 2L. The second half of the solution re- 
verses the twist produced by the first half. In all cases we have 
checked, solutions of this type are the lowest-energy solutions 
without net twist. 

A notable feature of Figure|2]and others like it is that the an- 
gles 9 and t/j quickly approach constant values away from the 
endpoints, which do not depend on the specific boundary con- 
ditions i22\ and i23\ . Thus we conclude that for long strips, 
the energy minimizing solutions have the following proper- 
ties: 9 and ip are constants such that 



sin 0; sin = A/L, 



(25) 
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in order to satisfy the global constraint. According to 
Eq. jl3a> . cos ip sin 9 — 0, which is compatible with Eq. ( 125 > 
only if 



■0 = 71/2, siii6' = A/L. 



(26) 



For solutions obeying Eq. ( I18> , the corresponding energy den- 
sity 



(27) 



is also constant, independent of s. Nevertheless, the solution 
oscillates because of Eq. (|8jl, which implies that 



(28) 



i?cos6' (A/L)2 
The period of oscillation is given by 



Sosc = 27ri?cos0 = 27ri?Vl - (X/L)^. (29) 

and the distance the solution travels when s traverses this arc- 
length is 

rA 



Zosc^2TTRsm9cos9 ^2nRy^l~(X/Ly-. (30) 

ij 

This last expression gives the period of oscillation one would 
measure upon taking a strip with L ^ R and pinning its ends 
at distance A in the laboratory. A three-dimensional visual- 
ization of such a strip appears in Figure|3jb). 

We conclude that long strips with metric gradients mini- 
mize their energy by curling up like telephone cords. They 
deviate from this quasi-uniform solution only when they must 
maneuver near the ends to obey boundary conditions, or form 
a kink in the middle to produce a solution with no net twist. 

Returning to the general case where Ci ^ C-z and writing 
out the Euler-Lagrange equations, we find that the uniform 
configuration i26\ is still a solution. The energy density is 
given hy £ = C2 tan^ 9/2R^ and all other features of the 
global minimum such as the period of oscillation remain un- 
changed. 



III. GEOMETRICAL CONSIDERATIONS 

The strip problem was originally motivated by experiments 
on two-dimensional sheets of deformed material. The ex- 
pectation was that by focusing upon a long thin strip, one 
could understand basic geometrical features of the full two- 
dimensional problem. Now that we have obtained a reason- 
ably detailed account of the ground state of the strip, it is clear 
that its significance for the original problem is limited. Long 
thin strips wrap up in a helix, like a telephone cord. This ge- 
ometry is impossible for the edges of a flower or a leaf. There- 
fore, in Section lrvl we return to the two-dimensional problem 
and search for some alternate approaches. As preparation, this 
section summarizes some basic facts concerning the differen- 
tial geometry of surfaces. 



Let r = [x, y, z) be a point in three-dimensional Euclidean 
space, where x, y, and z are the usual Cartesian coordinates 
with respect to a fixed laboratory frame. A surface may then 
be defined in parametric form by 



•{u, v) = (^x{u, v), y{u, v), z{u, v) 



(31) 



where u and v are parameters whose range is not specified 
for the moment. Throughout this section, derivatives will be 
abbreviated by di ~ d/du and 82 ~ d/dv. The elements of 
the metric tensor are then defined from 



gap = {dar) ■ {dpr). 



(32) 



A second fundamental form (tensor) is defined as follows: 
Let — fa (u, w) be the unit vector that is perpendicular to the 
surface at point (u, v) and varies continuously with u and v. 
The elements of the second fundamental form are then given 
by 



^3 • daC 



(33) 



The two symmetric tensors g^p and dap play a significant role 
in the theory of surfaces, as is apparent in standard textsQ 



A. Trumpets 

In order to gain familiarity with the types of surfaces stud- 
ied in Section Hvl we first consider an elementary example 
defined by 



f — p{v)[co&uei + sinue2) + C,{v)e^, 



(34) 



where ei, 62, and 63 are constant unit vectors along the three 
axes of the laboratory frame, while p{v) and C,{v) are func- 
tions of V alone, and are further restricted by the condition 



(35) 



where the primes indicate derivatives with respect to v. The 
surface is thus specified by the single function p = p{v). In 
particular, the metric tensor is then given by 



911= P , 922 = 1, 512 = 0. 



(36) 



In order to calculate the second fundamental form, we first 
construct the orthonormal trihedral 

fi = —dir — — sin uei + cos ue2 
P 

r2 — d2'r — p'(cosuei + sin 1*62) + \/l — p'^e^ (37) 
f3 = p'2(coswei + sinMe2) - p'e3, 

where fi and f 2 are tangent to the surface, while f 3 = fi x f 2 
is perpendicular The second fundamental form is then calcu- 
lated from Eq. (I33> to yield 



dll = -pVl - ^22 = 

dl2 = 6^21 = 



(38) 
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A special case of this class of surfaces is the ordinary cylin- 
der with unit radius, obtained with the choice p = 1 and 
( = V. The metric is then Euclidean (gu — 1 — 522,. 9i2 — 
921 — 0) and the elements of the second fundamental form 

are dn = -1, ^22 = 0, di2 = ^21 = 0. 



The corresponding orthonormal trihedral is then defined from 




FIG. 4: Rendering of the surface produced by the metric in Eq. <39> 
for It G [0, 27r],and v € {—00, 0] (in color). 



A more interesting example is obtained by the choice 

, 1 



l-v' 



(39) 



which leads to an axially symmetric surface with variable ra- 
dius p = p{v) and height z = ({v) calculated from Eq. j35t . 
This surface reduces to the ordinary cylinder in the limit 
V —00, but its radius increases with increasing v (or z). 
It would appear that the radius of this surface would eventu- 
ally grow to infinity as v 1. Actually, this limit cannot be 
obtained because of a singularity that develops earlier when 
p'^ = 1 or u = 0. The actual surface is illustrated in Fig- 
ure |4l we will refer to it as the trumpet. It starts as a cylin- 
der with radius p — 1 and terminates at a cusp with radius 
p ^ 2. Near the cusp the elements of the second fundamen- 
tal form calculated from Eq. ( I38> approach the characteristic 
limits dii 0, ^22 00, while di2 vanishes everywhere. 

This trumpet may be considered the most primitive model 
of a flower. However, a true flower is a trumpet that is allowed 
to grow beyond the cusp, and displays numerous ripples, be- 
cause biology and elasticity dictate slight modifications of the 
metric and second fundamental form needed to evade the sin- 
gularity. Geometry alone cannot resolve what happens as one 
tries to push a surface beyond such a cusp, but may help to 
classify the possibilities. 

Hence, in the remainder of this section, we formulate an 
inverse problem of sorts that will enable us to interpret the 
explicit results of Section lTvl 



B. Compatibility conditions 

The elementary example we have just considered is signifi- 
cantly generalized by considering surfaces for which the coor- 
dinate curves form an orthogonal net, which means that they 
are characterized by a metric of the form 



82? 



rs = ri X r2- 



'911 V-922 
We now consider the elementary integrability condition 



(41) 



did2r = d2dif 
^ c^i(\/522^2) = 92(^511^1) 

idly/g^)f2 + y/922dif2 = {d2y^)fi 



(42) 



Contracting both sides of Eq. ( I42t with the unit vectors fi, f2 
and fa in turn we obtain 



(fi • dif2) 



(f2 • d2fl) : 



d2^/9ii 
V922 
diy/g22 

Van 



(43a) 



(43b) 
(43c) 



We observe that for ^911 — I — v/R, 922 = 1, Eq. ( 143 a> 
reproduces Eq. As in Section [IIJ we parameterize the 
trihedral with the Euler angles of Eq. to obtain the three 
fundamental equations 



cos9di(j) + diijj - 



cos 6*920 + 92^' - 



d2^/gii 
di\/g22 



(44a) 
(44b) 



^(711 ( — sin 6 cos ijj d2(t> + sin i/i 820) 

= \/922(sin6'sin?/'<9i0 + cos?/'9i6'), (44c) 

which will provide the basis for subsequent development. It 
should be noted that the compatibility conditions Eqs. \AA\ are 
formulated entirely in terms of the metric tensor gap and that 
the elements of the second fundamental form do not appear 
explicitly. In fact, once a solution of Eqs. ( I44t is available, 
dap can be computed from Eq. ( I33t as 



dii = \/9ii{h ■ difi) 

= \/gii(^ — sin^cos V'Sic/) + smil)di9 

d22 = y/922{r3 ' 92^) 

— y/922 sin 9 sin -0 920 + cos tp 829 

dl2 = %/ff22(f^3 ■ dif2) 

= yigiJ (^sin0sin-0 9i0 + cos'ipdi9 
d2i = y/g^ii'Ts ■ d2fi) 

= V^ffn — sin 9 cos V' ^20 + sin i/" ^2 



(45a) 



(45b) 



(45c) 



(45d) 



511 = 9ll{u, v), 522 = 922{u, v), 512 = 0. 



(40) 



Note that the symmetry condition c?i2 = 1^21 is not explicit in 
Eq. ( I45t . but is enforced by Eq. ( I44c> . 



7 



As an elementary illustration, we return to the case of a 
trumpet characterized by a metric of the form ^/gu = p{v), 
JgTi — 1, and g\2 — 0. It is straightforward to verify that 



^ = 0, = 7r/2 + M, cos ( 



-p , sm( 



(46) 



is a solution of Eqs. ( I44> that reproduces the trumpet discussed 
earlier in this section. 

Our aim in Section BlI CI is to show that Eqs. J44t can be 
used as evolution equations actually to calculate trumpet- 
like solutions for the general class of metrics given by 
Eq. (|^) through the solution of an initial value problem. The 
connection with the standard Gauss-Codazzi compatibility 
conditions! 7] is briefly discussed in Section lTlIDI 



C. Evolution equations 

Consider the trumpet with metric given by Eq. ( I39t . When 
u ^ — oo, the surface approaches a cylinder. Can one choose 
some large negative value of v, suppose that the surface for 
this value of is a circle, and integrate towards v = know- 
ing nothing but the metric and reconstruct the trumpet? The 
answer is yes. One can find three first-order equations that ex- 
press the changes of Q, ip, and (j) with respect to v and integrate 
them forward as if u is a time variable. 

Finding the equations is not completely straightforward. 
Eqs. (I44> consist in three first-order partial differential equa- 
tions for the three Euler angles 0, ip, and (p. However since 
Eq. ( 144 a> involves derivatives in u only, algebraic manipula- 
tion alone does not allow one to solve for 820, 821^, and 82(j>- 

Note that derivatives of appear in Eqs. (I44> but not it- 
self. One can use the first two of Eqs. (I44> to express 8icj) and 
i920 in terms of the other two Euler angles. Removing (j> in 
this way is only permitted, however, if after solving for 8i(j) 
and one imposes the condition 8i82(t> = 828141 to obtain 



(V5ii .922 + 

where 



G 




g22Gcot9 



(47) 




One now solves the three Eqs. i44\ and Eq. ( I47> for 826, 

824^, 8i4> and 824- In expressing the results, it is convenient 
to use the expressions for the second fundamental form in 
Eq. ( 145 1 as shorthand for combinations of derivatives. One 
obtains 

di28i9 + Gcos'4 



826^ 



92-4 = -r- 




cos6di28i 



G sin -4 cot 6*] 
82-4 



cos 6, 



cos 9' 



(49a) 

(49b) 
(49c) 



These are the basic evolution equations for the Euler angles. 
Whenever 8i4 appears, it should be viewed as shorthand for 



di4) 



/g^8i4j + 82 



(49d) 



The procedure for constructing surfaces progresses as fol- 
lows: We specify the values of 6, ip, and </> for all u and some 
value of V — Vf). To construct a trumpet, we choose vq suffi- 
ciently negative so that we can use the expressions for cylin- 
der -0 = 0, 4 — tt/2 + u, and 9 = 7r/2 to set initial val- 
ues of the Euler angles for u G [0, 2tt]. Next, we calculate 
the right hand sides of Eqs. ( I49> using these initial values, 
and update each Euler angle through explicit Euler integra- 
tion, 9{v + dv) = 9{v) + 829 dv. It is easy to construct the 
surface f. To do so, form the trihedral from the Euler angles 
through Eq. (0. Then from Eq. ( I41> one has 
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(50) 



which means that if the surface is specified on the line v — vq, 
its future evolution can be determined as well. 

The integration process is extremely rapid. We have used 
it in order to reproduce the trumpet depicted in Figure |3 The 
process of integrating forward terminates at v = because 
dii vanishes there for all u, and two denominators in Eq. ( I49l l 
become singular 

We have also used the equations to generate surfaces from 
more general metrics, and a characteristic result is exhibited 
in Figure|5l Here we chose the metric 

522 = 1. (51) 

The surface has a slight three-fold modulation that is quite 
clear if one plots Euler angles, but rather faint when one looks 
at the surface itself. As in the case of the trumpet, the integra- 
tion process stops near to v — 0. Examining the reason, one 
finds that once again du vanishes. This three-fold trumpet 
does not have cylindrical symmetry, and du does not vanish 
simultaneously for all u. Instead, it first vanishes for three 
values of u. However, the evolution equations are singular 
and cannot be integrated past this point. Our original hope 
had been to use evolution equations to obtain flower-hke so- 
lutions such as the one displayed in Figure |8] Unfortunately, 
the rippling edge necessarily involves an oscillation in curva- 
ture, and du must oscillate in sign where ripples are visible. 
Thus our integration procedure is intrinsically unable to obtain 
solutions of this type. However, through trials with the evo- 
lution equations, we have found that they successfully create 
surfaces over ranges of u and v where du does not vanish. 

We note in closing that it is not difficult to generalize the 
equations of this section to the case where gi2 7^ 0. One takes 



{fil Jg22)82r ■ fl 



V'^- 912/911922 



(52) 



and otherwise proceeds as before. The resulting expressions 
are somewhat lengthy, and as we have not made use of them, 
we do not record them here. 
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symbolic form: 



FIG. 5: Three-fold trumpet produced by integrating Eqs. <49> for- 
ward from a cylindrically symmetric initial condition at vq = ~2n, 
using the metric in Eq. <51> (in color) . 



D. Gauss-Codazzi 



+ Ri 

dii 



92^11 = 51^12 + di^ {12} ^ ^^27 I } 

92^12 = did22 - d2j [^2} ^ '^^^ { 22 } ' 



(56) 



where { 2p } are the Christoffel symbols. Eqs. \56\ are again 
written in the form of evolution equations and can be solved 
with initial conditions starting at w ^ —oo given the ordinary 
cylinder values dn = —1, ^22 = and di2 = 0. Once a 
solution of Eqs. j56> is obtained at some finite v, the actual 
construction of the surface proceeds through the solution of a 
compatible system of linear equations: see Eqs. (39.8) ofRef. 
0. 



Our treatment has been based upon the compatibility con- 
ditions Eq. (I42> which are formulated purely in terms of the 
metric tensor gap, while the second fundamental tensor is a 
derived quantity displayed in Eq. ( I45> . By contrast, both ten- 
sors gap and dap appear as fundamental variables in the for- 
mulation of the Gauss-Codazzi compatibility conditions em- 
ployed in standard treatments IlZJ. It would then be interesting 
to demonstrate that the Gauss-Codazzi equations can actually 
be derived starting from the compatibiUty conditions Eq. J42> 
and the definitions Eq. ( I45> . 

As an example, we write the first two conditions in Eq. ( I42> 



ditp = — cos 9di(j) — 
d2ip = — cos 9d2(j) + 




(53) 



and impose the integrability condition 81821^ — d2diip to ob- 
tain 



sme{died2(t> - di(j)d2e) 
'diy/g^ 



di 




(54) 



On the other hand, a straightforward application of definitions 
(I45I together with Eq. dstt yields 



rfii'^22 — rfi2'^2i — -R1212, 

'di^/g^ 



^^1212 — —y/911922 



(55) 



dr-^)+d2r-^) 

^ Vdll ' ^ \/922 ' 



where R\2\2 is indeed the coiTect expression for the Rie- 
mann tensor element associated with a metric in the form of 
Eq. ( I40t . Thus Eq. ( 15 6> reproduces the celebrated Gauss equa- 
tion. 

We have not yet attempted a corresponding derivation of the 
two Codazzi equations starting from our Eqs. J42t and ( 145 > 
but merely state here all three Gauss-Codazzi equations in 



IV. GROUND STATE FOR FLOWERS 
A. Numerical technique 

We now return to techniques employed recentlyfs'] for the 
construction of surfaces through minimization of an elastic 
energy functional, and examine again the results in light of 
what we learned in Section HII] We begin by defining more 
carefully than has been done previously the problem that 
needs to be solved. In rough outline, one wants to take a thin 
flat sheet of material, impose a new metric gap upon it, and 
ask how it deforms in response. A more precise specification 
of the problem follows: 

Differential geometry describes a mapping between two 
spaces: a reference configuration described by variables 
(m, u), and a surface described by r(u,w). Experiments on 
deformed surfaces are performed by taking a flat sheet of ma- 
terial, deforming it in some controlled way, and then allow- 
ing the material to buckle in space. To construct a numerical 
model of the system we will first describe the reference state, 
coiTesponding to undeformed material, and then a discrete set 
of variables that describes a sheet of material moving about in 
three dimensions. 

The experimental reference state consists in a flat slab of 
material, much wider and longer than it is thick. Imagine 
therefore positions in a thin sheet described by (x, 2), where 
z g [0, t], and the thickness t is small. Pick iV points within 
this sheet, and label them by = (xi, zi), where i ranges 
from 1 to A^. In practice, we will take these points to sit on 
regular lattices, but they could be randomly distributed. Each 
point has a number of near neighbors: label these near 
neighbors with j e n{i). Describe the vector between two 
near neighbors by r^^ = — r^. We can now write down 
an energy functional which is constructed precisely so that 
its ground state gives back this reference configuration. This 
functional is defined on a new collection of variables fi, where 
again i ranges from 1 to A^, and the neighbor list j e n{i) is 
the same as in the reference configuration. Now, however, 
the points are free to move anywhere in three-dimensional 
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space. One can think of them as describing arbitrary defor- 
mations of the original thin sheet. If a particle i has particle 
j as a neighbor in the reference configuration, then particle j 
remains in the list of neighbors no matter how the sheet de- 
forms. 

The significance of neighbors is provided by an energy 
functional that depends upon the squared distance between 



pairs of neighbors. Define 



Uo 



4 



ri\ then 



(57) 



By construction this energy functional has the property that 
if every particle returns to the reference location f^, then 
the energy is zero. This ground state is not unique, for the 
energy is also unchanged if the locations of all the particles 
are rotated and translated in three-dimensional space, as when 
one picks up a piece of cardboard and translates and rotates it. 
Depending upon details involving the numbers of neighbors 
of each particle, there may be additional degeneracies in the 
ground state as well, but we will not worry about this point 
right now. 

Our numerical model for deforming the sheet is to go to 
each bond and stretch it so that the square distance between 
near neighbors i and j adopts the new value Ify Thus we have 
the energy functional 



U = 



hi) I 



(58) 



and direct minimization of U has been employed to obtain 
most of the three-dimensional figures in this paper. 

In order to connect Eq. ( I58t with the discussion in the pre- 
vious sections, we must explain the connection between bond 
lengths lij and metrics. We use the following prescription. 
For each bond ij, choose values of u and v through 



(59) 



In other words, take u and v to be the x and y coordinates 
of the midpoints of bonds in the reference configuration. We 
define a target metric through three functions 



(60) 



5ii(w,w),522(">«):and gl^iu^v). 



One could choose, for example the functions in Eq. J36> if 
one wanted to recover a trumpet. Using the target metric, one 
obtains a new equilibrium length squared tf^ for the bond be- 
tween points i and j through 



Q/3 



^13 



(61) 



So, for example, if a bond in the reference configuration lies 
along the x axis at position (u, v) given by Eq. \59\ , its new 
length is g\^{u, v). The reason that we call g^^ the target 
metric rather than the metric is that one can put any collection 
of particle locations one wishes into the functional Eq. j58L 
not just particle locations that correspond to surfaces with 
metric g^^. Thus the metric of the surface obtained through 



numerical minimization may in principle be different from the 
target. We will be attempting to determine whether the ground 
states of U produce surfaces whose metric equals the target 
metric. 

One might wonder why U involves squares of bond lengths. 



rather than (r^ — Uj) . The answer is that Eq. ( 15 8> leads to 
conventional nonlinear elasticity in the continuum limit, while 
the alternative does not. To obtain the continuum limit, let 
r(w, v) be a continuous vector field, and write 

nj ~ (rf, • V)r. (62) 



Recalling Eq. ( I32> . substitute Eq. \62\ into Eq. ( I58t to obtain 



E-Oa/ t \-0l3 

rij [gap - ga(})r,- 



Q/3 



(63) 



that is, the energy is given by subtracting the target metric 
from the actual metric, and vanishes when the two are equal. 
From Eq. ^ one sees that the appropriate generalization of 
the Lagrangean strain tensor to situations with target metrics 

is 



(64) 



If the target metric is a unit tensor, E reduces to the con- 
ventional Lagrangean strain tensor of nonlinear elasticity, and 
when deformations are small it reduces further to the strain 
tensor of linear elasticity. One can write 



E 

a/3 



ij 



-Op 

•r- 
13 



(65) 



For a particular arrangement of mass points, one can perform 
the sums over and obtain a specific quadratic functional 
depending upon the components of E^. We do not need 
these expressions here are will not pursue them further 

However, we will spell out the particular reference config- 
uration that has been used to produce results for this paper. 
It consists in either one or two layers of a triangular lattice. To 
be completely explicit, use three integers Imn to describe the 
point locations rather than the single index i: 



r«„„ = (l,0,0)? + (i,^,0)m+(^,0, 



(66) 



where / and m range over positive and negative integers, and n 
equals to produce a single layered structure, or ranges over 
and 1 to produce a two-layered structure. In the two-layered 
structure, each particle has nine nearest neighbors, each at unit 
distance; six with the same value of n on the same horizontal 
sheet, and three with a different value of n on a different hori- 
zontal sheet. In these crystalline structures, the near neighbors 
of particle i are all the particles j for which jr^^ P = 1. 

We record one final technical point about the numerical 
techniques. Given particle locations r^, we will want to view 
the particles as describing a continuous surface, and to con- 
struct its metric. To do so, focus on the lower sheet of particles 
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FIG. 6: Enumeration of neighbors surrounding point i, used to de- 
scribe the construction of numerical metrics in Eq. <68> . The diagram 
shows the particles sitting in the reference configuration fi described 
by Eq. |66}. 



(n = in Eq. (I66» and label the neighbors of particle i as in 
Figure|6l Denote by Sr^ the square distance from point i to the 
neighbor located at in the figure and so on for the remaining 
five neighbors. Then we have numerical representations of the 
metric 



511 



522 



512 



Svq + Sri 
2 

Sri + Srj + Sri + Sr?^ - Srl/2 - Srl/2 



(67) 



Sri ^ '^^2 ^ '^^4 + '^^5 



2V3 



B. Numerical experiments 

The questions we wish to pose about target metrics are: 

1 . What are the ground states of Eq. (|63}? 

2. When do these ground states correspond to smooth sur- 
faces in the continuum limit? 

3. When the ground state is a smooth surface, under what 
conditions does the target metric equal the metric of the 
surface? 

Observe that ground states of Eq. ( I63> always exist. The 
functional is positive definite, and for a finite number of par- 
ticles must have one or more global minima. If the reference 
configuration does not have bending stiffness, however, the 
resulting ground state corresponds to a nondifferentiable sur- 
face. To illustrate this point, we use Eq. ( I39> as a target metric, 
and work on the domain 



u e [0,27r], and v e [-37r,0]. 



(68) 



We represent the system with a reference crystal 200 columns 
long, 346 rows high, but only 1 layer thick (in Eq. ( I66t . 



FIG. 7: Image of surface created by minimizing Eq. <63t for a refer- 
ence crystal with 200 x 346 x 1 particles, and target metric given by 
Eq. <69t . The surface achieves the target; the total energy summed 
over all particles is less than 10^^. However, the solution is nowhere 
smooth, and resembles tree bark. It can be understood as an embed- 
ding of the target metric whose first derivative is nowhere continuous 
(in color). 



I e [0,200],m e [0, 346],n = 0). Since the reference con- 
figuration is infinitely thin, there is no source of bending stiff- 
ness. Using the target metric 



522 = 1 



1 



(69) 



we minimize U, and the result is displayed in Figure Q The 
total energy U has converged below 10~^, and each bond has 
reached its target value to better than two parts in 10''. How- 
ever, the surface is not smooth, nor is the configuration dis- 
played in the figure plausibly unique. As we know from Sec- 



tion llll Al there do exist smooth surfaces whose metric equals 
this target metric, but the minimization routine under these 
conditions does not find them. Venkataramani et al. fT3l have 
pointed out the existence of non-smooth surfaces of the sort 
appearing in the figure. 

Low energy states of Eq. j63t look very different when 
one moves from reference crystals with one layer to reference 
crystals with two layers (In Eq. i66\ . n E [0, 1]). Now the 
system possesses some stiffness, and pays an energy penalty 
for bending too rapidly. Once again minimizing Eq. ( I63> for 
the parameter range in Eq. ( I68t . we find the minimum energy 
state that is essentially the trumpet shown in Figure|3 The en- 
ergy of this structure is 74.3; it cannot converge to zero, since 
springy elastic material has been wrapped into a cylinder. 

We compared the metric g as computed through Eq. ( I68t 
with the target metric for the structure in Figure |4] For 
each metric component, the difference is featureless, and ap- 
proximately equal to— SxlO^'^at every point, due to a slight 
compression needed to bend the inner layer of the reference 
crystal into a cylinder We conclude from the numerical cal- 
culations that the metric and target metric are equal within 
numerical accuracy. We emphasize that in Section Ull CI we 
reproduced the trumpet in Fig. [^through evolution equations; 
the metric and boundary conditions suffice to determine the 
surface in this case. 

As a more interesting exercise in direct minimization of 
Eq. i63\ . we extend the domain beyond the critical point and 
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FIG. 8: Image of surface created by minimizing Eq. <63> for a ref- 
erence crystal witli 200 x 346 x 2 particles, and target metric given 
by Eq. <69t . corresponding to a domain where v varies from — Stt to 



Once V > 0, trumpet solutions given by Eq. J34> no longer exist, 
and the evolution equations in Eq. <49> are incapable of finding so- 
lutions. The minimum energy state is the smooth flower-like surface 
with seven-fold symmetry displayed here. As shown in Figure|9| the 
metric of this surface does not equal the target metric. Creation of 
this surface required long series of minimizations. The process be- 
gan by placing particles in a cylinder with a Euclidean metric, and 
very slowly changing the metric until it reached the desired target 
value, continually minimizing the functional Eq. <63> along the way. 
Attempts to find the surface more quickly resulted in higher-energy 
structures with creases (in color). 



look for a solution in 



ue[0, 27r], and uG[-37r, i]. 



(70) 



The significance of increasing the range of v is that the the- 
ory of Section 1111 Al is unable to find a smooth surface for 
V > 0. In contrast, upon minimizing U once more, we find the 
flower-like surface shown in Figure|8] Now when we subtract 
the target metric Eq. ( I69> from the metric actually achieved, 
Eq. ( I68t . the difference is visible, as shown in Figure |9] The 
seven-fold pattern in the surface is reflected in seven-fold os- 
cillations in gii and 522. The off-diagonal component of the 
tensor, gi2, remains zero within numerical accuracy. We con- 
clude that in this case, no smooth surface is able to reproduce 
the target metric, and numerical minimization finds a met- 
ric close to the target that is capable of producing a smooth 
surface even for v > 0. We see no reason that this surface 
should be considered unique. In all likelihood, its details de- 
pend upon the thickness of the sheet. According to Audoly 
and Boudaoud|4], one should expect the surface to be increas- 
ingly ramified as its thickness diminishes. 

As a final exercise, we took the metric components com- 
puted numerically and depicted in Fig. |5J inserted them into 
the evolution equations Eqs. i49i and attempted to reproduce 
the structure in Figure|8l This attempt was only partially suc- 
cessful. The evolution equations were unable to proceed past 
points where dn approached zero. It is clearly possible for 
dii to vanish without a surface developing singular cusps at 
that point, but our numerical routines would require very del- 
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FIG. 9: Metric component deviations for the flower depicted in Fig- 
ure[8| Each panel shows a contour plot of the difference between the 
metric component gap and the target metric g^^. The off-diagonal 
metric element gi2 is not illustrated because its numerical values are 
too small to be discerned at the scale of the figure (in color). 



icate cancellations in order to proceed past such a point. 



V. CONCLUSIONS 

Given a surface, it is completely straightforward to compute 
the metric. The inverse problem of finding surfaces compati- 
ble with a given metric is much more difficult. In setting out 
on the studies recorded in this paper, we had two goals. First, 
we wanted finally to determine the low-energy configurations 
of long strips with linear gradients in metric. This we have ac- 
complished. Second, we hoped to determine conditions on the 
relatively simple metrics thought to create the shapes of leaves 
and flowersi 1] that would enable surfaces to be reconstructed. 
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In the second task we have been only partially success- 
ful. We found evolution equations enabling construction of 
surfaces from initial conditions and metric alone. However, 
at points when dn, a component of the second fundamental 
form, vanishes, the equations become singular and cannot be 
integrated further. We are unable to tell when this singularity 
really reflects the impossibility of creating a surface compat- 
ible with the metric, and when it simply reflects a technical 
defect in the method of construction. 

Some intuitive understanding of this situation can be ob- 
tained by thinking about a piece of paper, with flat Euclidean 
metric. Imagine holding the bottom of the paper completely 
straight, along the u axis, while the left edge of the paper runs 
along the v axis. There is an infinite number of shapes the 
paper can take, which correspond to all different bends possi- 
ble at different points along v around axes parallel to u. All 
of these solutions have vanishing curvature in the u direction; 
dii =0, and the evolution equations i49i accordingly are un- 
able to make any predictions. That is, in some cases the van- 
ishing of dii can correspond to a genuine uncertainty, based 
upon initial conditions, concerning how the surface should 



evolve. 

On the other hand, for the surface in Figure|SJ the curvature 
dii oscillates in sign and therefore passes through zero, yet the 
surface exists. Therefore, the vanishing of du can be compat- 
ible with the existence of a surface even if our purely geo- 
metrical methods cannot find it beyond points where rfn = 0. 
Minimization of Eq. ( I58> still produces a surface because elas- 
ticity resolves the questions that cannot be answered by geom- 
etry. 
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